Spin dynamics of hole doped Y 2 BaNiOs 
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Starting from a multiband Hamiltonian containing the relevant Ni and orbitals, we derive an 
effective Hamiltonian H e ff for the low energy physics of doped Y2BaNiOs. For hole doping, H £ ff 
describes O fermions interacting with S=l Ni spins in a chain, and cannot be further reduced 
to a simple one-band model. Using numerical techniques, we obtain a dynamical spin structure 
factor with weight inside the Haldane gap. The nature of these low-energy excitations is identi- 
fied and the emerging physical picture is consistent with most of the experimental information in 
Y 2 _,Ca,BaNi0 5 . 

PACS. 75.10J - Heisenberg and other quantized localized spin models. 

PACS. 75.25 - Spin arrangements in magnetically ordered materials (neutron studies, etc). 
PACS. 75.50L - Spin glasses. 



The properties of antiferromagnetic (AFM) chains 
have been a major topic in the study of low dimen- 
sional magnetism. Haldane's conjecture about the exis- 
tence of a gau.in spin 5=1 one dimensional (ID) Heisen- 
berg systemsEl has inspired a considerable amount of 
researchHa. Recently, there has been wide interest in 
Y 2 BaNi0 5 , a near ideal realization of the 5=1 AFM 
Heisenberg chaup, and the effects of different dopants 
on specific heatQ and magnetic propertiesotl. While re- 
placement of Ni 2+ (5=1) by non-magnetic Zn 2+ or Mg +2 
simply cuts the spin chains, the replacement of off-chain 
Y 3+ by Ca 2+ introduces holes primarily on oxygen sitesQ. 
In the last case, resistivity decreases by several orders of 
magnitude, and inelastic neutron scattering reveals new 
states inside the Haldane gapQ. In addition, the magnetic 
behavior of the Ca doped system, markedly different from 
the Mg-doped one, is sapm-glass like with a characteristic 
temperature T g ~ 2K.B 

It has been shown that spin models with localized 
site or bond impurities lead to in-gap states for cer- 
tain parameters&El. While the results are interesting, 
these models can only be valid for completely local- 
ized added holes, which seems not to be the case in 
Y2-j;Ca«BaNi05l3. Models for mobile holes have been 
studiedU'Ej, but with effective hole hopping t more than 
one order of magnitude smaller than realistic valuesliil 
In Ref£3, the essential physics is reduced to that of a 
weakly m obil e 5=1/2 spin (representing a "Zhang-Rice 
" doublett3ta) with exchange interaction J' = with its 
nearest-neighbor (NN) Ni spins. In the static limit, this 
model is known to lead to-p-bound state inside the Hal- 
dane gap when J' is weaUJii3. However J' results strong 
and ferromagnetic (FM)c3 and, in agreement with Ref.t3, 
we find no bound states in this static case. Thus, the 
origin of the observed states inside the Haldane gap in 
Y2- x Ca :E BaNi05L3 is not yet clear. 

In this Letter, we start from a multiband Hamilto- 
nian H m b for MO5 linear systems, containing all essen- 



tial orbitals and interactions. An effective low-energy 
Hamiltonian H e f f is derived. By exact diagonalization of 
H e ff, we provide an explanation of the observed inelastic 
neutron spectrumQ. In addition, our results are consis- 
tent with the measured gap, x-ray absorption spectrumQ, 
magnitude of the intra- and inter-chain exchanged, mag- 
netic susceptibility and muon spin relaxatiora. 

The starting Hamiltonian is the extension to,-the ID 
MO5 system of the one used by van Elp et alEi for a 
NiC>6 cluster, restricting the basis of Ni orbitals to the 
r 2 and 3d x 2_ y 2 (z is the direction along the chain): 
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where cj Q(T creates a hole with spin a on the orbital a of 
site i. The interactions included are the on-site Ni repul- 
sion Ud, Ni Hund interaction J H (for O sites = 0), O 
on-site repulsion U p and Ni-0 NN repulsion U p d, for the 
six O atoms surrounding a Ni one. For simplicity, the 
basis of O orbitals considered is reduced to the 2p Zl and 
the two linear combinations of 2p orbitals lying in the xy 
plane having optimum hybridization with their NN Ni or- 
bitals (the amount of other O 2p states, in the low energy 
manifold, is of the order of 1%X From atomic data we 
obtain J H = 1.6eV. From RcfEJ U d = U( 2 E g ) = lOeV. 
The valus&iof U p = 4eV and U p d — 1.2eV were taken 
from Ref.li-i The hopping parameters were taken from 
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Ref.B, scaling the Ni-0 (0-0) hoppings with an r -7 / 2 
(r~ 2 ) law, and multiplying them by a common factor /. 
Finally, / and A, the energy necessary to take a hole 
from Ni 2+ and put it in any of the six nearest 2 ~ ions, 
were chosen to giwe J = 0.023eV and a charge transfer 
gap E g — 1.84eVEil in agreement with experiments!! 

The resulting values / = 0.85 (implying a 3d^ z 2_ r 2 — 
2p z hopping of 1.61eV ) and A = 7.4eV are reasonable 
taking into account that A should increase to the left of 
the periodic table, and A = e p — + U p d ~ 4.8eV in the 
cupratesEJ. To gain further confidence on the resulting 
parameters, we have compared other two quantities with 
experiment. Solving exactly i? m h, including the repul- 
sion Vd (V p ) between a Ni 2p 3 / 2 core hole and a Ni 3d 
(NN O 2p) hole in a Ni06 cluster, we have calculated the 
shift towards lower energy in the Ni L3 x-ray absorption 
spectrum in going from polarization along .the chain to 
a perpendicular one. Taking Vd — V^=10eVt-lwe obtain 
0.80eV, while the experimental shift is O.91eV0. We have 
also calculated the exchange along b direction (perpen- 
dicular to the chains) Jt,, starting from the exact solution 
of two NiOg clusters, and treating the hoppings connects 
ing the two clusters with the cell-perturbation methodE-3. 
We obtain J& = 0.173AT = 6.4 x 10j~ 4 J, near the upper 
bound 5 x 10~ 4 J estimated in Ref.tl ._. _ 

The low-energy reduction procedurelliliijO which we 
found most convenient fox. the present problem is the one 
used successfully in RefE3 to reproduce low-energy pho- 
toemission and magnetic properties of the three-band 
model for the cuprates. The form of the low-energy 
effective Hamiltonian H e ff is determined by a canoni- 
cal transformation. For hole doping H e ff has the form 
of a spin-fermion H s f (or Kondo-Heisenberg) model as 
in Ref.EJ. For electron doping, H e ff has the form of 
a one-band model Hu (Eq. (3) with J' = and 
t = — 0.351eV). The procedure eliminates linear terms in 
the Ni-0 hopping, retaining second order contributions 
and the fourth-order exchange J. However, the values of 
the parameters are determined fitting the energy levels 
of H s f to the corresponding ones of H m b in conveniently 
chosen clusters (N12O11 for electron-doped and NiOg for 
hole-doped systems). H s f can be written as: 
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Here i labels a Ni site, and i + S denotes its two NN O 
atoms along the chain, p\ + g a creates an effective 2p z hole 
(it contains information of other Ni and O orbitals) of 
spin a at site i + S, and Si is an effective spin 

1 (1/2) at site i (i + 5). We obtain t x = 0.64eV, t 2 = 
0.89eV and J K = 1.40eV; e = 8.97eV is irrelevaet for 
the spin dynamics, but determines the charge gapcil. 
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FIG. 1. (a) Full (dashed) line: dynamic structure factor 
S(u!,tv) (5 , (lj,7t)/5 )of 0(NiOs)io system with open bound- 
ary conditions and one (zero) added hole, (b) Correlation 
functions E(ui), M(ui) and L(u>) defined in the text. 

In Fig. 1(a) we show the dynamical structure factor 
S(q,u) = i / e iuJt dt (S 9 (t).S_ g (0)> at wave vector q = w, 
in a (NiOs)ioO system with one added hole described 
by H s f, using open boundary conditions (OBC), and 
compare it with the result for the undoped system. At 
first sight, OBC appears to be an inconvenient choice for 
studying states inside the Haldane gap, since an open 
chain has always low-energy spin excitations associated 
with the nearly free spins at the ends&Qtl However, the 
use of periodic BC leads to spurious low-energy excita- 
tions due to the frustration generated by effective FM 
interactions introduced by the hole dynamics. This is 
particularly clear when the spin-spin correlation between 
the ends of the open chain is FM. Hence, we used OBC 
and identified the excitations related to the end states 
by the dynamical response of the system to a staggered 
magnetic field decaying exponentially along the chain. 
For no added holes, S(tt, uj) shows a peak at energy u> ~ 
J = 23meV corresponding to the Haldane gap-, (shifted 
from u> ~ 0.4 J due to the finite size of the chaina), and a 
peak at to ~ 3meV corresponding to excitations localized 
at the ends of the chain. This is the behavior expected 
for Zn or Mg doping. S(-k,oj) for one added hole (rep- 
resenting the Ca doped system) is markedly different. 
The Haldane peak shifts to lower energy, its intensity de- 
creases, and significant spectral weight appears near the 
middle of the Haldane-jgap. Both features are in agree- 
ment with experimentQ. Among the three peaks below 
the Haldane gap, the second one lies nearly at the same 
energy as that of the low energy peak of the undoped 
system. One would then suspect that it is mainly related 
with the ends of the chain and it should disappear in the 



2 



thermodynamic limit, since Ca doping does not break the 
chains. This suspicion is confirmed by inspection of the 
correlation function E(u) = f dt e luJt (§ \O(t)O(0)\ \) 
where O = Y^j e ~ 1+ and \m) is the component 
with S z — m of the 5=3/2 ground state. E(u>) measures 
the response of the system to an excitation localized at 
one end of the chain. The result is shown in Fig. 1(b). 
There is an intense structure around 3meV consisting of 
three peaks. The first and the third one are practically 
absent in 5(tt, co) for symmetry reasons. The second cor- 
responds precisely to the second peak of S(ir,ui). 

What is the origin of the other two peaks below the 
Haldane gap? Guided by the wave function of the 
ground state and the first few excited states, we have 
looked for dynamical correlation functions with large 
intensity at these peaks to unravel their nature. The 
peak at 9meV (near the middle of the Haldane gap) is 
related to a smaller effective magnetic interaction be- 
tween the Ni spins at a distance of one lattice param- 
eter and a half from the O hole and their NN Ni spins. 
To check this and since the O hole tends to be lo- 
calized near the middle of the chain, we have calcu- 
lated M{lo) = Jdt e iwt (§|O 2 (t)O 2 (0)|§) with 2 = 

Su_ c i>3/ 2 (— I) 5 e%~\i- c \Sj, where the half-integer c de- 
notes the center of the chain. As shown in Fig. 1(b), 
M(ui) has only an intense peak at 9meV, and a smaller 
contribution at the Haldane peak. Replacing 2 by 
, the corresponding dynamical correla- 
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tion function L(w) has an intense peak only at the en- 
ergy (around 1.7meV) of the lowest peak in S{ir,uS), and 
smaller contributions at the other low-energy peaks (see 
Fig.l (b)). Note instead that the responses to the end ex- 
citation E(uS) and to the excitation peaked at the middle 
of the Haldane gap M(ui) have negligible contributions 
at 1.7meV. The energy of the lowest possible excitation 
depends non-monotonically on system size (for example 
it lies at 0.6meV for 0(NiOs)8), but we find always an 
excitation at very low energies. These results suggest the 
presence of weak magnetic links near the O hole and the 
presence of a FM polaron which is a bounded state be- 
tween the charge and a local magnetic excitation. The 
origin of this polaron, can be described as follows: the 
large values of t%, t 2 , and Jk (relative to J) tend to 
stabilize a doublet between an effective Ni spin at any 
site i and a holo4n-t|he orbital pi+s + Pi-s (" Zhang- Rice" 
doublet (ZRD))ll30, and induces a residuaLEM interac- 
tion between the ZRD and its NN Ni spinsEi This FM 
interaction can be interpreted in terms of a NN attractive 
interaction between the ZRD and a magnetic excitation. 
This attractive interaction gives place to a bound state 
(small polaron) with a binding energy which is a frac- 
tion of J. Within this picture, the effect of the operator 
M(ui) is to make a spin excitation inside the polaron, 
while the effect of L(u>) is to make an excitation in the 
spin which is beside the polaron. It is clear that the AFM 
interaction between the latter spin and the one inside the 
polaron must be reduced, giving place to a weak link, by 



the effects of t±, t 2 , and Jk which tend to polarize the 
spin background. This explains the origin of the very 
low-energy excited state, generated by the application of 
L(u>) to the ground state. 

We .discuss briefly a further reduction of H s f to a 
simpleEl one band Hamiltonian Hit, since. in..general, 
it is successful in the case of the cupratestrBB. The 
low-energy physics is dominated by "Zhang-Rice" dou- 
blets. The overlap of some of these many body states 
is too large (1/3 compared with 1/8 in the case of the 
cupratesE3) rendering cumbersome-to obtain Hu as an 
expansion in powers of the overlaps We then construct 
orthogonal Wannier functions -Ki at each Ni site, write 
H sf in this basis using p i+s = (2/tt) I]„(-l)' i+1 ( 2n - 
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calculate P\H s fPi (P\ is a projector over 



locaLorthogonal Zhang- Rice doublets), and map it into 
H\$3. .-Unfortunately, contrary to the case of the 
cupratestlj, P\ does not commute even approximately 
with H a f in any limit. Local quadruplets are impor- 
tant in the low energy physicsEa. Assuming for simplicity 
Jk = t\ + t 2 , and including only the two largest terms 
(besides J) we obtain: 



t Pj,]+2S ( Sj.Sj+2S + ^ ) - J' Sj.Sj+2S 



(3) 



+ J J]] Sj.Sj+i 



Here j runs over the positions of the Zhang- Rice 5=1/2 
spins Sj and Pjj±2S permutes Sj with its nearest 5=1 
spins Sj±28- The resulting values of t and J 1 are t = 
0.36ti + 1.027* 2 ~ leV - 40J, and J' = ~ 
0.25eV, in agreement with sirriple expectationsE3llj. The 
S(tt,uj) obtained using Hu, is very different from that of 
H s f and displays only one peak (around 6meV) below 
the Haldane gap. This result is in fact more similar to 
the response of the undoped system. 
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FIG. 2. (a) Top: schematic representation of the ground 
state of H e f for one added O hole (indicated by the small 
up arrow). The dashed line indicates a region in which the 
interchain interactions Jd are frustrated due to the change of 
the interchain short-range magnetic correlations produced by 
the O hole. (b)The same for the first excited state. The Ni 
states with S z = 1, 0, —1 are indicated by j", 0, 1 respectively. 
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To shed light onto the underlying physics, it is instruc- 
tive to solve H s f in a M3O2 chain. The ground state has 
S=5/2 and is nearly degenerate to an 5=3/2 state. Ex- 
cept for the parity of the wave functions, the result can 
be interpreted as if the central static ZRD were strongly 
FM coupled to one of the neighboring Ni spins and very 
weakly FM coupled to the other one. This quasidegener- 
ancy is not obtained with Hu for any set of parameters, 
indicating that the inclusion of quadruplets is crucial to 
have this weak link. This image is confirmed by the anal- 
ysis of the dominant terms in the wave function of the 
ground and first excited states represented in Fig. 2, 
and provides a plausible explanation of the susceptibility 
measurements of Kojima et al.u: at temperatures above 
the effective exchange of the weak link, the chains behave 
as being cut, and a segment between two ZRD's can have 
zero, one or both ZRD's FM coupled with it, with prob- 
abilities 1/4, 1/2 and 1/4 respectively. If the number of 
Ni atoms in the segment is even, the spin of the segment 
is zero. Then, the corresponding total spin of the seg- 
ment plus coupled ZRD's is 0, 1/2, for the three cases 
respectively. If the number of sites in the segment is odd, 
the spin of the segment is equal to 1 and the correspond- 
ing total spin is 1, 3/2 and 2 respectively. This gives an 
average (S(S + 1)) = 17/8 in excellent agreement with 
the observed Curie constants. 

An important feature of the ground (first excited) state 
of H s f in the 10-site chain, is that the end spins display 
FM (AFM) correlations. When the interchain exchange 
Jf, is included, this implies that for 10% Ca doping, near 
five Ni spins have the wrong sign of the magnetic cor- 
relations with neighboring chains, implying a frustration 
energy of roughly 20 Jb = 3.2K (see Fig. 2). Since this 
energy is withiruthe order of magnitude of the randomly 
(due to disorder!] ) distributed lowest excitation energies, 
a spin glass behavior with characteristic temperatures 
near 3K is expected, in agreement with experimental. 

In summary, starting from a multiband Hamiltonian 
similar to a one used previously for NiOEJ but including 
only 3d e g orbitals, we have derived a low-energy Hamil- 
tonian for doped Y2BaNi05, and used it to understand 
magnetic properties of the Ca doped system, particularly 
the additional states inside the Haldane gap. Contrary 
to previous predictionsEJ, we obtain that these states are 
quite different to these produced by Zn doping and sub- 
stantial weight near the middle of the Haldane gap is 
obtained only for Ca doping. 
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